AD-A266  010 


TECHNICAL  REPORT  ARCCB-TR-930I6 


AD 


GENERALIZED  MANDELBROT  RULE 
FOR  FRACTAL  SECTIONS 


L.V.  MEISEL 


APRIL  1993 


US  ARMY  ARMAMENT  RESEARCH, 
DEVELOPMENT  AND  ENGINEERING  CENTER 

CLOSE  COMBAT  ARMAMENTS  CENTER 
BENET  LABORATORIES 
WATERVLIET,  N.Y.  12189-4050 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


93-13951 

I  (fill  tfiii  »>•  ..... .  ^  ■ 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  authorized 
documents . 

The  use  of  trade  naae(s)  and/or  manufacturer(s)  does  not  constitute 
an  official  indorsement  ot  approval. 


DESTRUCTION  NOTICE 

For  classified  documents,  follow  the  procedures  in  DoD  5200. 22-M, 
Industrial  Security  Manual,  Section  11*19  or  DoD  5200. 1-R,  Information 
Security  Program  Regulation,  Chapter  IX. 

For  unclassified,  limited  documents,  destroy  by  any  method  that  will 
prevent  disclosure  of  contents  or  reconstruction  of  the  document. 

For  unclassified,  unlimited  documents,  destroy  when  the  report  is 


no  longer  needed.  Do  not  return  it  to  the  originator. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-Q 1 88 


rdoeiting  burden  tor  thu  collection  of  ntotmiiion  estimated  to  average  '  nour  per  tewaon-ve.  including  the  time  lot  reviewing  invtructiom.  veefcnmg  ,.s!  -*'.i  source* 

gathering  «nd  maintaining  the  data  needed,  and  competing  and  reviewing  the  collection  of  information  Seng  comments  reg+tamg  ttu*  burden  estimate  or  $nv  nxer  jsoect  of  ?n«s 
collection  of  information,  inclining  lu^emom  for  reducing  this  burden  to  Washington  needauaner*  Services.  Directorate  for  information  Operations  400  fteoons.  i 4  S  ,*tfer*>n 
Oavrs  Highway,  Suite  U04.  Arlington,  va  22202*4302.  and  to  the  Office  of  Management  and  0udget.  Paperworn  Reduction  Project  (Q 704-0 1 98).  Washington,  -JC  20SC3 


1.  AGENCY  USE  ONLY  (Leave  blank)  I  2.  REPORT  DATE 


3.  REPORT  TYPE  AND  OATES  COVERED 


April  1993  Final 


4.  TITLE  AND  SU8TITLE 

GENERALIZED  MANDELBROT  RULE  FOR  FRACTAL  SECTIONS 


S.  FUNDING  NUMBERS 

AMCMS:  6 16224 H 1800 
PRON:  1A12186EA11A 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  ADDRESS(ES) 
U.S.  Army  ARDEC 
Benet  Laboratories.  SMCAR-CCB-TL 
Watervliet,  NY  12189*4050 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARCCB-TR-93016 


9.  SPONSORING/ MONITORING  AGENCY  NAME(S)  ANO  ADORESS(ES) 

U.S.  Anny  ARDEC 

Close  Combat  Armaments  Center 

Picatinny  Arsenal,  NJ  07806-5000 


10.  SPONSORING  i  MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

Published  in  Physical  Review  AI5  (American  Physical  Society) 


IZa.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited 


13.  ABSTRACT  (Maximum  200  words) 

Mandelbrot's  rule  for  sections  is  generalized  to  apply  to  the  Hentschei  and  Procaccia  fractal  dimension  at  arbitrary  q  3nd  on 
arbitrary  sections.  It  is  shown  that  for  almost  all  (n-m)-dimensional  sections,  D°(q)  =  D°‘“(q)  +  rn,  where  the  D'(q)  are  box¬ 
counting  Hentschei  and  Procaccia  generalized  fractal  dimensions  of  r-dimensional  sections  of  homogeneous  fractal  point  sets 
in  R°.  The  rule  is  shown  to  apply  for  finite  "thickness"  sections  as  well  as  ’true"  sections.  A  more  general  form  of  the  rule 
applicable  to  inhomogeneous  fractal  sets  is  also  presented. 


1 4Trac^s.TcIeneral ized  Dimensions.  Box-Counting,  Fractal  Dimensions 


15.  NUMBER  OF  PAGE5 

8 


16.  PRICE  COOE 


17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 


N5N  7S40-01 -280-5500 


18,  SECURITY  CLASSIFICATION  I  19.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 


OF  THIS  PAGE 

UNCLASSIFIED 


OF  ABSTRACT 

UNCLASSIFIED 


Standard  form  298  Pev  2-89) 

O'etcnoed  by  'itd  <39**8 

298-102 


TABLE  OF  CONTENTS 


ACKNOWLEDGEMENTS .  it 

INTRODUCTION .  1 

THEORY  .  1 

SUMMARY  AND  DISCUSSION  .  4 

REFERENCES .  5 


List  of  Illustrations 

1.  Illustration  of  the  box  and  strip  definitions  in  R: .  6 


lion  f-pr 

|  NUS  CRA 41 

i  one  tab 
1  Unannounced 
Justification 


J 


U 

□ 


By _ 

Distribution  I 


A va liability  Codes 


Dist 


Av^il  j-)g/#r 

Special 


ACKNOWLEDGEMENTS 


The  author  gratefully  acknowledges  P.J.  Cote  and  Mark  Johnson  of  Benet  Laboratories  for 
stimulating  discussions  and  encouragement  concerning  this  work. 


a 


INTRODUCTION 


Mandelbrot’s  rule  (ref  1)  states  that  for  almost  all  (n-1) -dimensional  sections 

DhB  -  »HB  ♦  1 

where  the  Drar  are  Hausdorff-Bouligand  dimensions  of  r-dimensionai  sections  of  homogeneous  fractal  sets 
in  R°.  This  important  rule  comprises  the  basis  for  a  substantial  body  of  numerical  fractal  analysis. 

The  following  generalization  of  Mandelbrot’s  rule  is  demonstrated  here:  For  almost  all 
(n-m)-dimensional  sections 


D\q)  =  Dn  m(q)  -  m  (1 

where  the  Dr(q)  are  box-counting  Hentschel  and  Procaccia  (refs  2,3)  generalized  fractal  dimensions  of 
r-dimensional  sections  of  homogeneous  fractal  point  sets  in  R“.  Equation  (1)  obtains  for  arbitrary  q  and 
n  >  m  5  1.  The  q  =  0,  m  =  l  form  of  Eq.  (1)  is  the  box-counting  formulation  of  the  standard 
Mandelbrot  rule. 

THEORY 

Let  the  point  set  in  R“  be  represented  as  points  in  an  arbitrarily-oriented  Cartesian  coordinate 
system.  Define  an  ensemble  of  (n-m)-dimensionai  ’’sections," 

TJJ)  *  r„«WI<  •  !•*>;  «> 


s  \x\X.fJ)  <  ij  s  XfJ)  +  6  ,  i  =  l  jn\  -t  <  x}  ^  L\  j  ~  m  +  \ji)  (2) 


where  L  is  large  enough  to  cover  a  substantial  portion  of  the  fractal  set,  5  is  "small,"  and  the  (Xj(J)}  are 
such  that  the  Tnjn(J)  fill  a  region  of  Rn  without  overlaps  or  gaps.  We  refer  to  Tnjn(J)  sections  as  ’’J-strips." 

Define  sets  of  boxes  (hypercubes)  of  side  E  aligned  along  the  coordinate  axes  in  Rn.  Let  each 
hypercube  in  R"  be  identified  by  the  position  of  one  of  its  corners  in  units  of  E,  i.e.,  the  box  at 
(I(l),...,I(n))E  is  indexed  [1(1), ...,I(n)].  Such  sets  of  boxes  define  corresponding  sets  of  boxes  of  side  E  and 
"thickness"  5“  in  each  J-strip:  the  box  corresponding  to  the  [1(1),..., I(n)]  box,  which  intersects  the  J-strip, 
induces  the  box,  which  is  indexed  [J;I(m+  l ),..., I(n)j,  in  the  J-strip.  There  will  be  approximately  E/5 
J-strips,  which  intersect  a  given  box  in  Rn.  Figure  1  illustrates  the  definitions  in  R2. 

Box-counting  fractal  analysis  is  generally  applied  to  large,  but  finite,  subsets  of  fractal  sets.  The 
fractal  subsets  are  assumed  to  be  large  enough  to  yield  converged  box-counting  dimensions.  Small  5 
implies  that  5  <  <  E^,,  where  is  the  smallest  hypercube  edge.  The  present  results  also  apply  for 
complete  fractal  sets  and  in  the  limit  5  -*  0. 

Let  Pi(i),„j(n)  denote  the  occupation  probability  for  the  [I(l),...,I(n)[  box  of  side  E  in  Rn.  which  is 
defined  as  the  number  of  members  of  the  subset  of  the  fractal  point  set  in  the  [I(l),...,I(n)[  box  divided  by 
the  number  of  members  of  the  subset  of  the  fractal  point  set  in  all  the  boxes  in  R". 

Let  N,J(m>l)  J(n)  be  the  number  of  points  of  the  fractal  subset  in  the  [J;I(m+l),...,I(n)|box. 
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Let  Sj  be  the  number  of  members  of  the  subset  of  the  fractal  point  set  in  the  J-strip: 

Sj  =  T1  NJAm+ Dr-M  ‘ 

Sj  is  proportional  to  the  J-strip  probability  cr,  in  Ra,  where 

Oj  s  Sj/N 


an^  ~  m  Sj  ~  E  l),... 


An) 


is  the  total  number  of  members  in  the  fractal  subset.  The  occupation  probability  for  the  [J:I(m+  l),...,I(n)J 
box  is  s  Njj(m+I)_j(n/Sj.  By  definition 

52  7lJAm~V),...An)  =  1  ■ 

Hm+l),...An) 


The  partition  function  Zn(q,E)  in  R°  is  defined  as 


Z'*(^v£)  a  E  Jin)  ■ 

m....An) 


(3a  i 


Since  the  J-strip  box  corresponding  to  the  [I(l),...,I(n)j  box  in  R°  has  volume  Q,  =  En(5/E)m, 


Furthermore,  each  box  in  R”  is  counted  approximately  (E/5)ra  times  in  the  sum  over  strips.  Thus,  the 
partition  function  Zn(q,E)  in  R“  can  be  expanded  as 


Z\qJE)  =  £  £ 


6 

m 

E 

£ 

6 

\* 


£  a)  Zn-m(J\q£) 


(3c) 


where  Znra(J;q,E)  is  the  J-strip  partition  function.  The  approximations  become  exact  in  the  limit  E  -*  0 
with  E/5  >  >  1  and  are  "good"  approximations  in  the  spirit  of  box-counting  algorithms  for  "small  enough" 
E  with  E/5  >>  1;  exact  equalities  obtain  in  the  E  -*  0  limit  for  true  sections  also. 
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Box-counting  generalized  fractal  dimensions  are  defined  as 


.  Wj£*g»l 

1  d!n(E) 


(4a) 


Dn'm(J;q)  . 

1  81n(£)  J,. 


(4b) 


where  the  right-hand  sides  are  usually  determined  by  least  squares  fitting  over  a  range  of  "small"  yardstick 
values  centered  about  a  yardstick  E<,. 

Since  ln(Z°(q,E))  and  ln(Z°'m(J;q,E))  scale  like  (q-l)D°(q)ln(E)  and  (q-l)Do  m(J;q)ln(E), 
respectively,  Eq.  (4)  implies  that 

Dn(q)  =  <Dn'm(q)>£g  +  m  <5) 


where 

<D’(q)>.  =  £  Joj  /  '£ZV'-,qJE^,  . 

J  J1 


(6) 


E0, 


(Since  the  strip  partition  functions  Zntn(J;q,E0)  are  of  the  same  order  of  magnitude  for  given  q  and 


</>'(«)>  =  £  DV,q)o,J  I  £  o], 

J  j' 


(6’) 


might  be  a  useful  approximation.) 

N.B.,  Eqs.  (5),  (6),  and  (6')  may  be  interpreted  for  inhomogeneous  fractal  sets. 

Excluding  "exceptional  sections,"  the  Dr(J;q)  are  independent  of  J  for  homogeneous  fractal  point 
sets  and  Eq.  (5)  is  equivalent  to  the  generalized  form  of  Mandelbrot’s  rule  for  box-counting  fractal 
dimensions  at  arbitrary  q  and  almost  all  (n-m)-dimensional  J-strips 

Dn(q)  =  Dn~m(J;q)  +  m  ,  for  almost  all  J  .  (7) 

The  existence  of  "exceptional  sections"  related  to  spec.al  symmetries  in  fractal  point  sets,  as 
discussed  by  Mandelbrot  (ref  1),  necessitates  the  inclusion  of  the  "almost  all"  caveat  in  the  statement  of 
the  rule  for  sections. 

The  steps  in  the  derivation  of  Eqs.  (5)  through  (7)  hold  in  the  8  -*  0  limit.  Thus,  the  principal 
results  hold  not  only  for  finite  "thickness"  (finite  5)  sections,  but  they  also  hold  for  "true”  sections  of  fractal 
sets. 
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SUMMARY  AND  DISCUSSION 


Mandelbrot’s  rule  relates  tbe  fractal  dimension  (for  q  =  0)  of  a  fractal  point  set  in  R°,  D°(0),  to 
the  fractal  dimension  (for  q  =  0)  of  subsets  of  the  fractal  point  set  in  (n-l)-dimensional  sections,  Dn  ‘(0). 
The  principal  results  of  the  present  work,  contained  in  Eqs.  (5)  through  (7),  generalize  Mandelbrot’s  rule 
for  sections  to  apply  for  box-counting  Hentschel  and  Procaccia  fractal  dimensions  D'(q)  at  arbitrary  q  for 
r-dimensional  sections  of  small  "thickness"  5  with  n  >  r  >  0.  Small  8  implies  that  S  <  <  E„,  where  E0  is 
defined  after  Eq.  (4).  Furthermore,  the  equations  hold  in  the  5  -»  0  limit,  and  thus  applv  for  'true’’ 
sections. 

One  might  work  with  finite  5  sections  to  analyze  point  sets  generated  by  nonlinear  mappings  or 
the  distribution  of  galaxies.  One  might  work  with  "true"  two-dimensional  sections  of  three-dimensional 
surfaces  in  a  multifractal  generalization  of  the  slit-island  method  (ref  4)  (which  determines  D(0)  of 
fracture  surfaces  by  perimeter-area  analysis  of  fracture  surface  sections). 

Since  sparse  regions  of  a  fractal  subset  dominate  the  average  in  Eq.  (6)  for  q  <  0,  numerical 
problems  in  box-counting  evaluations  of  D°(q)  are  likely  at  q  <  0  Modified  box-counting  strategies,  which 
put  less  stress  on  sparsely  occupied  strips,  might  improve  fractal  analysis  at  negative  q. 
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[1(1),  1(2)]  box 


[J;  1(2)]  box 


Figure  l.  Illustration  of  the  box  and  strip  definitions  in  R‘.  The  box  indexed 
(If  1),I(2)|  of  side  E,  whose  ’volume"  is  E:  in  R\  induces  the  box 
indexed  [J;I(2’)J  of  side  E  in  the  J-strip.  The  J-strip  is  a 
one-dimensional  section  of  "thickness"  5. 
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